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Abstract — A feedback scheme for preparation of photon 
number states in a microwave cavity is proposed. Quantum 
Non Demolition (QND) measurement of the cavity field provides 
information on its actual state. The control consists in injecting 
into the cavity mode a microwave pulse adjusted to increase the 
population of the desired target photon number. In the ideal 
case (perfect cavity and measures), we present the feedback 
scheme and its detailed convergence proof through stochastic 
Lyapunov techniques based on super-martingales and other 
probabilistic arguments. Quantum Monte-Carlo simulations 
performed with experimental parameters illustrate convergence 
and robustness of such feedback scheme. 

I. Introduction 

In [10], [5], [4] QND measures are exploited to detect 
and/or produce highly non-classical states of light trapped in 
a super-conducting cavity (see [6, chapter 5] for a description 
of such QED systems and [1] for detailed physical models 
with QND measures of light using atoms). For such experi- 
mental setups, we detail and analyze here a feedback scheme 
that stabilize the cavity field towards any photon-number 
states (Fock states). Such states are strongly non-classical 
since their photon numbers are perfectly defined. The control 
corresponds to a coherent light-pulse injected inside the 
cavity between atom passages. The overall structure of the 
proposed feedback scheme is inspired of [3] using a quantum 
adaptation of the observer/controller structure widely used 
for classical systems (see, e.g., [7, chapter 4]). The observer 
part of the proposed feedback scheme consists in a discrete- 
time quantum filter, based on the observed detector clicks, to 
estimate the quantum-state of the cavity field. This estimated 
state is then used in a state-feedback based on Lyapunov 
design, the controller part. In theorems [T] and |2] we prove the 
convergence almost surely of the closed-loop system towards 
the goal Fock-state in absence of modeling imperfections 
and measurement errors. Simulations illustrate this results 
and show that performance of the closed-loop system are 
not dramatically changed by false detections for 10% of the 
detector clicks. In [2] similar feedback schemes are also 
addressed with modified quantum filters in order to take 
into account additional physical effects and experimental 
imperfections. [2] focuses on physics and includes extensive 
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closed-loop simulations whereas here we are interested by 
mathematical aspects and convergence proofs. 

In section HIl we describe very briefly the physical system 
and its quantum Monte-Carlo model. In section Hill the feed- 
back is designed using Lyapunov techniques. Its convergence 
is proved in theorem [T] Section |IV] introduces a quantum 
filter to estimate the cavity state necessary for the feedback: 
convergence of the closed-loop system (quantum filter and 
feedback based on the estimate cavity state) is proved in 
theorem|2]assuming perfect model and detection. This section 
ends with Theorem [3] proving a contraction property of the 
quantum filter dynamics. Section[V]is devoted to closed-loop 
simulations where measurement imperfections are introduced 
for testing robustness. 

The authors thank Michel Brune, Serge Haroche and Jean- 
Pierre Raimond for useful discussions and advices. 

II. The physical system and its jump dynamics 




Fig. 1. The microwave cavity QED setup with its feedback scheme (in 
green). 

As illustrated by figure [H the system consists in C a 
high-Q microwave cavity, in i? a box producing Rydberg 
atoms, in i?i and i?2 two low-Q Ramsey cavities, in D an 
atom detector and in S* a microwave source. The dynamics 
model is discrete in time and relies on quantum Monte-Carlo 
trajectories (see [6, chapter 4]). It takes into account the back- 
action of the measure. It is adapted from [5] where we have 
just added the control effect. 

Each time-step indexed by the integer k corresponds to 
atom number k coming from B, submitted then to a first 
Ramsey 7r/2-pulse in Ri, crossing the cavity C and being 
entangled with it, submitted to a second 7r/2-pulse in R2 
and finally being measured in D. The state of the cavity 
is described by the density operator p^. Here the passage 
from the time step fc to fc + 1 corresponds to the projective 



measurement of the cavity state, by detecting the state of the 
Rydberg atom number k after leaving R2. During this same 
step, an appropriate coherent pulse (the control) is injected 
into C. Denoting, as usual, by a the photon annihilation 
operator and by = a^a the photon number operator, the 
density matrix pk+i is related to pk through the following 

D{ak)MkPkMlD{-ak) 



jump-relationships: pk+i 



Tr( MfcPfcM+) 



where 



the measurement operator Afj. — Mg (resp. Afj. — M^), 
when the atom k is detected in the state \g) (resp. |e)) 
with 



Me 



(1) 



Such measurement process corresponds to an off- 
resonant interaction between atom k and cavity where 
(f)R is the direction of the second Ramsey 7r/2-pulse (i?2 
in figure [TJ and is the de-phasing-angle per photon. 

• The probability Pg ^ (resp. Pe,k) of detecting the atom 
k in \g) (resp. |e)) is given by Tr[MgPkMg) (resp. 

Tr(MePfcMe). 

• D{ak) is the displacement operator describing the co- 
herent pulse injection, D{ak) = cxp{ak{a^ — a)), and 
the scalar control ak is a real parameter that can be 
adjusted at each time step k. 

The time evolution of the step fc to A: + 1, in fact, consists 
of two types of evolutions: a projective measurement and 
a coherent injection. For simplicity sakes, we will use 
the notation of p^^i, to illustrate this intermediate step. 
Therefore, 



Pfc+i = D(afc)pj,^iZ)(-Qfe) (2) 



In the sequel, we consider the finite dimensional approx- 
imation defined by a maximum of photon number, In 
the truncated Fock basis (|"))o<,i<„raax, N corresponds to 
the diagonal matiix (diag(?i))Q<jj<jj,™x, p is a {rf" + 1) x 
(n°"" + 1) symmetric positive matrix with unit trace, and 
the annihilation operator a is an upper-triagular matrix with 
(v^)i<,i<n'™ upper diagonal, the remaining elements 
being Q. We restrict to real quantities since the phase of any 
Fock state is arbitrary. We set it here to 0. 

III. Feedback scheme and convergence proof 

We aim to stabilize the Fock state with n photons char- 
acterized by the density operator p = |n) To this end 
we choose the coherent feedback ak such that the value of 
the Lyapunov function V{p) = 1 — Tr (pp) decreases when 
passing from p^+k to pk+i- Note that, for a small enough, 
the Baker-Campbell-Hausdorff formula yields the following 
approximation 



D{a)pD{—a) ^ p — a[p, a' — a] + ^ [[p, a' ~ a],a^ — a] (3) 



up to third order terms. Therefore, for ak small enough, we 
have 



Tr(^D(afc)pj^^iD(-afc)P^ 

Tr (pfc+ 1 p) - QfeTr ([p^^i , a+ - a]p) 

+ ^Tr ([[p^^i, at- a], at -a]p 

Thus the feedback 

ak = ciTr (^[p,a^ - a]pk+^'j 
with a gain ci > small enough ensures that 

Tr(ppfc+i) - Tr (pp^+i) > ^Iti' ([p,a'' - a]p^^i) I , (5) 



(4) 



since Tr([pfe+i,at 



a\p 



-Tr 



thermore, the conditional expectation of Tr \^pPk^ 
ing Pk is given by 



Fur- 
know- 



EfTrfp, 



PPk+i 



since [p, Mg] 
Thus 



+ P.,fcTr 



pMgPkMl 

Pg.k 

\ Pe,k 



and AfJAfg 



= Tr (ppfc) 
MlM, = 1. 



E(Tr(ppfe+i) I pfc)>E(Tr 



PPk- 



Pk = Tr (ppfe) 



and consequently, the expectation value of V{pk) decreases 
at each sampling time: 



EiVipk+i))<E{Vipk)) 



(6) 



Considering the Markov process pk, we have therefore shown 
that V{pk) is a super-martingale bounded from below by 0. 
When V reaches its minimum 0, the Markov process pk 
has converged to p. However, one can easily see that this 
super-martingale has also the possibility to converge towards 
other attractors, for instance other Fock states which are all 
the stationary points of the closed-loop Markov process but 
with V{p) — 1 instead of 0. Following [9], we suggest the 
following modification of the feedback scheme: 



ciTr (jp, at 



argmax Ji ipD{a)pf,^iD{~a 



ifVipk) <l~£ 
if V(pfc) > 1-e 
(7) 

with ci,e,a > constants. 

Theorem 1: Consider (|2|i and assume that for all n G 
{0, . . . , n™"} we have + 7^ mod {tt/2) and that 



# i cos 



+ n(t>] \ ne {0, 



} = + 1. 



bR + i 



Take the switching feedback scheme ^ with a > 0. For 
small enough ci > and e > 0, the trajectories of (|2]| 
converge almost surely towards the target Fock state p. 

Remark 1: The second part of the feedback (|7]l, dealing 
with states near the bad attractors, is not explicit and may 
seem hard to compute. Note that, this form has been partic- 
ularly chosen to simplify the proof of the Theorem [T] and 



in practice, one can take it to be any constant control field 
exciting the system around these bad attractors and ensuring 
a fast return to the inner set. 

Remark 2: The controller gain ci can be tuned 
in order to maximize at each sampling time fc, 

Tv(^D{ak)Pk+^D{-ak)pj for pf,^i neai- p. Up to 
third order term in Pk^i — p, © yields to 



Tr 



+ 



(Ti([p,a'^ -a]p^^^ij^ (^Ci-^Ti(^[p,a'' -a][p,a'' -a]j^ . 

Thus ci = 1/Tr([p,at - a][p,at - a]) « l/(4n + 2) for 
^ n implies a maximum decrease at the sampling time, 
up to third-order terms in pk — p. 

In order to prove the Theorem [T] we need some classical 
tools from stochastic processes namely the Doob's inequality 
and the Kushner's asymptotic invariance Theorem [8]. These 
results are been recalled in the Appendix. 

Proof of Theorem [7] It is divided in 3 steps: in a first 
step, we show that for small enough e and by applying the 
second part of the feedback scheme, the trajectories starting 
within the set {p | > 1 — e} reach in one step the set 

{p I V{p) < 1 — 2e} and this with probability 1; next, we 
show that trajectories starting within the set {p \ V{p) < 
1 — 2e}, will never hit the set {p \ V{p) > 1 — e} with a 
uniformly non-zero probability p > 0; finally, we will show 
that, the trajectories of the quantum filter converge towards p 
for almost all trajectories that never hit the set {p \ V{p) > 
1 — e}. This is then an immediate conclusion of the Markov 
property that the trajectories of the quantum filter with the 
feedback scheme (|7]i will converge almost surely towards p. 

Step 1 : We start by considering the process starting on the 
level set {p \ V{p) = 1}. We have the following lemma: 

Lemma 1: Consider p a well-defined density matrix such 
that Tr {pj5) = 0. We have 



mm max 



Tr (pD(a)M,pMj"D(-a)) 



Tr M, 



spMl^ 



> 0. 



We denote any argument of the above min-max problem by 

a{p) e [-a,a\. 



Proof of Lemma [7} Define 



s £ {g,e}. 



Tr(M3pA/J) ' 

The matrices Mg and Me being diagonal, we trivially have 
Tr (psp) = 0. Let us fix s and assume that for all a G 
[-«,«], 

Tr{pDia)p,D{-a))^0. (8) 

Decomposing ps as a sum of projectors we have ps = 
J2"kLi ^k^s Hk,s) {i^k,s\ , where Afe^^ are strictly positive 
eigenvalues and ip^^s are the associated normalized eigen- 
states of Ps {m — 1 corresponds to the case where ps is a 
projector). The equation ([8]l, clearly, implies 

{i>k,s \ D{~a)n) ^ 0, Vfc G {I,--- ,m},Va G [-a,Q]. (9) 

Fixing one k G {1, • • • , m} and taking = V'fe.s, noting 
that D{—a) = exp(— a(a^ ~«)) ™d deriving j times versus 
a around we get 



(V' I (a^-a)%) =0, Vj^O,...,/!" 



But the family ((a^ — '^)"'")o<j ^^^^ rank. This is a 
direct consequence of [1 1, Theorem 4]. It is proved there that 
the truncated harmonic oscillator ^ \(j>)^ = —{iN+v{t){a^ — 
a)) , is completely controllable with the single scalar 
control v{t). If the rank r of ((at — a)P <„m" 
strictly less that n'"" + 1, then according to Cayley-&amilton 
Theorem the rank of the infinite family ((a^ — a)^ l"-))p>o 
is also r. Take |^), a state orthogonal to this family. For 
any control v{t), the state starting from \n) remains 
orthogonal to |^). Thus it will be impossible to find a control 
v{t) steering \<j))^ from \n) to |^). 
Since the rank of ((at — a)^ \n)) 



0<p<n"' 



IS maximum. 



( [Tol l implies \4>k) = and leads to a contradiction. □ 

Applying the compactness of the space of density matri- 
ces, we directly have the following corollary: 
Corollary 1: There exists an e > such that 



lr{pD{a{p))MspMlD{-a{p))) 
mi — ^ -. ^ > 2e 

pe{Tr(pp)<E} jj. / 



spMt^ 



(11) 



for s — g,e and where a{p) is defined in Lemma [T] 
Proof of Corollary [7} We take 



5 = inf min 

pG{Tr(pp) = 0}sS{9,e} 



Tr {pD{a{p))MspMtD{-a{p))) 



Tr M, 



spMl'j 



By Lemma [T] and the compactness of the set {p | Tr (pp) = 
0}, we know that 5 > 0. By continuity of Tr (pp) with respect 
to p and by compactness of the space of density matrices, 
there exists 7 > such that 

. Tr {pD{a{p))MspM!D{-a{p))) S 
mi mm — y r > — . 

pe{Tr(pp)<7}aS{s,e} TiIMsPM}^ 



Therefore, by taking e = min(7, (5/4), clearly, ( fTTI ) holds 
true. □ 

Through this corollary, we have shown that whenever the 
Markov process hits the set {Tr (pp) < e}, it is immediately 
rebounded to the set {Tr (pp) > 2e} and this with probabil- 
ity L 

Step 2: Let us assume that the process starts within the 
set {Tr (pp) > 2e}. 

Lemma 2: Initializing the Markov process within the set 
{p I V{p) < 1 - 2e}, Pk will never hit the set {p | V{p) > 
1 — e} with a probability p > > 0. 

Proof of Lemma ]2} By (|6]l, the process V{pk) is, clearly, a 
supermartingale. One only needs to use the Doobs inequality 
(cf. Appendix, Theorem |4]l and we have 



P{ sup V{pk) >l-e)< 

0<k<oo 



Vipo) 



2e 



(10) 



and thus p > 1 - (1 - 2e)/(l - e) = e/(l - e). □ 

We have shown that starting within the inner set 
{Tr (pp) > 2e} there is a uniform non-zero probability 
of e/(l — e) for the process, to never hit the outer set 
{Tr(pp)<e}. 
Step 3: 

Lemma 3: The sample paths pk remaining into the set 
{Tr (pp) > e} converge in probability to p as fc ^ 00. 



Proof of Lemma ^ Consider the function W(/9) = 1 

\2 



Tr [pp) . For s 



and similarly 

>V(pe) 



= 1 



^ Tr{pMl-pM,f 
Tr(M,pMt)" 

Tr (MgpMj)' 
Tr (MepMl)' 



Tr(pp) = 



-Tr(pp)^ 



(12) 



(13) 



Furthermore, whenever a is given by the first part of the 
feedback law, we have 

W{D{a)pD{~a)) - \^{p) < -2eci |Tr (^[p, a"!" - a]p) |^ (14) 

where we have applied Q together with the fact that 

\Tx{D{a)pD{-a)p)\ + \Tx{pp)\>2e 

since p is inside the set {^^{pp) > e}. Apply- 
ing ([2]), (fTSI) . ( fTsl l and (fT4l i for the paths never leaving the 
set {Tr {pp) > e}, we have 

E(>V(pfc+i) I Pk)-W(Pk) < -2eci|Tr([p,at -a]p^_^i) |' 



1 COS ^ 




1* 

1 


sin 




1* 


Tr 




Tr 


{^MePkMp^ 





1 Tr (pfep)2 



Noting thatXr (AfgPfcM]) > q, Tr (MePfcAft) > 0, 

Tr {MgPkMl) + Tr (M^PkM^) = 1 and by Cauchy- 
Schwartz inequality, we have 



1 cos I ij3i±^ + flC 



IxiMgPkMl 



/ I cos f + nc, 



Tr [A'IgPkMl 



TriUePkMl 



+ 



sin I + 



Tr ( M^pkMl 



(Tr ( MgPfcMj ) + Tr ( M^pfcA/l ) ) > 



+ n(f> ) + sin 



+ nq 



= 1, 



with equality if and only if Tr (A/gp^Aft j 



2 / it>a+4' 



COS' I + n0 1 . We apply, now, the Kushner's invari- 

ance Theorem (cf. Appendix, Theorem |5]l to the Markov 
process p^ with the Lyapunov function W{pk)- The process 
Pk converges in probability to the largest invariant set in- 
cluded in 

{p I Tr (M.pMl) = cos^ (^^^ + ,1^) } 

fl {p I Tr ([p, at - a]Af,pAft) = 0, s = 5, e} . 
In particular, by invariance, p belonging to this limit set 
implies Tr(A/,pAft) = ^^/(^J^t^"^ for s ^ g,e. 



Taking s = g, and noting that Alg = Afj, this leads to 
Tt{M^p) = Tr{M'^pf . However, by Cauchy-Schwartz 
inequality, and applying the fact that p is a positive matrix, 
we have Tr (Af^p) = Tr {M^p) Tr (p) > Tr (Af ^p)^ , with 
equality if and only if M^p and p are co-linear. Since Mg has 
a non degenerate spectrum, p is necessarily a projector over 
one of the eigen-state of Af^, i.e., a Fock state \n), for some 
n G {0, . . . , Finally, as we have restricted ourselves to 
the paths never leaving the set {p | Tr (pp) > e}, the only 
possibility for the invariant set is the isolated point {p}. □ 
Lemma 4: pk converges to p for almost all paths remain- 
ing in the set {Tr (pp) > e}. 

Proof of Lemma ^ Define the event P>£ = 
{ijj e I pfe never leaves the set {Tr (pp) > 

e}}. Through Lemma [3] we have shown that 

limfc^ooP(||pfc-p|| >^ I P>.) = 0, > 0. By 

continuity of V{p) = 1 — Tr (pp), this also implies that 
limfc^ooP(1^(pfe) > 5 I P>,) = 0, V(5 > 0. As V{p) < 1, 
we have 

E{V{pk) I P>.)<f{V{pk)>5\ P>,) 

+ Sil-FiV{pk)>S\ P>,)). 

Thus limsup;,^^E(V"(pfc) | P>e) < (5, V(5 > 0, and 
so limk^ac>E{V{pk) I P>e) = 0. By Theorem El we 
know that V{pk) converges almost surely and therefore, 
as V is bounded, by dominated convergence, we obtain 

E(iimk^^vipk) I p>e) = o.n 

Now, we have all the elements to finish the proof of the 
Theorem[T] From Steps 1 and 2 and the Markov property, one 
deduces that for almost all paths pk, there exists a K such 
that Pk for k > K never leaves the set {Tr (pp) > e}. This 
together with the step 3 finishes the proof of the Theorem. □ 

IV. Quantum filtering for state estimation 

The feedback law (|7]i requires the knowledge of Pk+i- 
When the measurement process is fully efficient and the jump 
model ^ admits no error, it actually represents a natural 
choice for quantum filter to estimate the value of p by p°" 
satisfying 



pt+i = D{ak)pt+iDi-ak) 



Tr(^M,,p-Ml) 



(15) 



where Sk = g ot e, depending on measure outcome k and 
on the control Uk- 

Before passing to the parametric robustness of the feed- 
back scheme, let us discuss the robustness with respect 
to the choice of the initial state for the filter equation 
when we replace pj^^i by p^'^i in the feedback (|7]i. Note 
that. Theorem [T] shows that whenever the filter equation is 
initialized at the same state as the one which the physical 
system is prepared initially, the feedback law ensures the 
stabilization of the target state. The next theorem shows that 
as soon as the quantum filter is initialized at any arbitrary 
fully mixed initial state (not necessarily the same as the initial 



state of the physical system (O) and whenever the feedback 
scheme d?) is appUed on the system, the state of the physical 
system will converge almost surely to the desired Fock state. 

Theorem 2: Assume that the quantum filter ( fTSl l is initial- 
ized at a full-rank matrix pg and that the feedback scheme (|7| 
is applied to the physical system. The trajectories of the 
system (|2|, will then converge almost surely to the target 
Fock state p. 

Proof of Theorem |2} The initial state pg' being full-rank, 
there exists a 7 > such that — 7^9 + (1 — 7)Po, where 
Po is the initial state of (|2| at which the physical system 
is initially prepared and pg is a well-defined density matrix. 
Indeed, p'" being positive and full-rank, for a small enough 
7, (pg — 7pg)/(l — 7) remains non-negative, Hermitian and 
of unit trace. 

Assume that, we prepare the initial state of another identi- 
cal physical system as follows: we generate a random number 
r in the interval (0, 1); if r < 7 we prepare the system in 
the state pg and otherwise we prepare it at pg. Applying 
our quantum filter ( fTSI l (initialized at pg') and the associated 
feedback scheme, almost all trajectories of this physical 
system converge to the Fock state p. In particular, almost 
all trajectories that were initialized at the state pg converge 
to p. This finishes the proof of the theorem. □ 

The quantum filter ( fTsT i admits also some contraction 
properties confirming its robustness to experimental errors 
as shown by simulations of figures [3] and |4] where detection 
errors are introduced. We just provide here a first interesting 
inequality that will be used in future developments. 

Theorem 3: Consider the process (|2]l and the associated 
filter ( fTSl ) for any arbitrary control input (ak)'kLi- We have 
E (Tr (pfcp-)) < E (Tr {pk+iPt+i)) , V/c. 
Proof Before anything, note that the coherent part of the 
evolution leaves the value of Tr(pfep5!') unchanged: 

Tr (pfe+iPfc'+i) = Tr (^D{ak)pi^^ipl\iD{-ak)^ 

= Tr(p,+ ,p-;,). 

Concerning the projective part of the dynamics, we have 



EfTr 



Pfc+iPfc+i j I Pfc,Pfc 



^ Tv{MsPkMtMsptMl) 
2^ 7 7^ ■ (16) 



Tr 



(m^p-'mJ) 



Applying a Cauchy-Schwarz inequality as well as the 
identity M^Mg + M^M^ = 1, we have 



E 



Tr (Ml M.pfcMjM.pg') 
Tr (^MsPtMl^ 



Tr (m,p-'MJ) 
J2 ^Tr {ulM^pkMlMsp'^ j (17) 

— g,e / 



> Txipkpt)- (18) 



Applying ( fT6l ) and ( [TtI i. we only need to show that 

i J Tr (m1 Ms pk Ml Ms pf^ 

\s — g,e / 

Noting, once again, that M^Mg+M^Me ~ 1, we can write: 
Tr(pfePfe')=E E Tr(M]M.PfeMtM,pl'). (19) 

s—g,e r—g.e 

and therefore ( fTSI ) is equivalent to 



^Ti(^MsPkMlMspfMl^ J2 \jTr{MrPkMlMrPtMl^ 

s — g,e r — g,e 

>Y1 E Ti-(A4M,pfeMtM.prV (20) 



rPk 

s — g,e r — g,e 

Note that as p^ and pf! are positive Hermitian matrices, their 
square roots, y/pk and y^pf^, are well-defined. Once again 
by Cauchy-Schwarz inequality, we have 

Tr [hdlMsPkMlMrPk^ = Tr [y^MlMs^^MlMr^'^ 
< jTiiv^MlMs^^MlMsy^' 



Tr ( \^m1 Mr Vp^^M 1 Mr y^' 

= ^Tr (^MsPkMlMsPtMl^ ^Tr (^MrPkM:^MrPtM:l.^ . 

Summing over s,r E {g, e}, we obtain the inequality (l20b 
and therefore we finish the proof of the Theorem [5] □ 

V. Monte-Carlo simulations 

Figure (|2|l corresponds to a closed-loop simulation with 
a goal Fock state n — 3 and a Hilbert space limited to 
n""" — 15 photons, pg and pg' are initialized at the same 
state, the coherent state exp(v^(a^ ~a)) |0) of mean photon 
number n. The number of iteration steps is fixed to 100. The 
dephasing per photon is (p — j^. The Ramsey phase is 
fixed to the mid-fringe setting, i.e. -\- n(f> = j. The 

feedback parameter with p^^^i instead of Pk+i) are as 
follows: ci — ^ 



4n+l' 



and a 



j_ 

10' 



Any real experimental setup includes imperfection and 
error. To test the robustness of the feedback scheme, a false 
detection probability 77/ = is introduced. In case of false 
detection at step k, the atom is detected in g (resp. e) whereas 
it collapses effectively in e (resp. g). This means that in (fTsl l. 
Sfc — g (resp. Sk — e), whereas in dU, it is the converse 
Mfe = Me (resp. Mk = Mg). Simulations of figure |3] differ 
from those of figure |2] by only rjf = j^: we observe for 
this sample trajectory a longer convergence time. A much 
more significative impact of rjf is given by ensemble average. 
Figure |4] presents ensemble averages corresponding to the 
third sub-plot of figures |2] and [3] For rjf —0 (left plot), we 
observe an average fidelity Tr (pfcp) converging to 100%: 
it exceeds 90% after fc = 40 steps. For 77/ = 1/10, the 
asymptotic fidelity remains under 80% and reaches 70% after 
30 iteration. The performance are reduced but not changed 
dramatically. The proposed feedback scheme appears to be 
robust to such experimental errors. 



VI. Conclusion 
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Fig. 2. A single closed-loop quantum trajectory in the ideal case (n = 3). 



In [2] more realistic simulations are reported. They include 
nonlinear shift per photon {N(f) replaced by a non linear 
function ^{N) in ([T]i) and additional experimental errors 
such as detector efficiency and delays. These simulations 
confirm the robustness of the feedback scheme, robustness 
that needs to be understood in a more theoretical way. In 
particular, it seems that the quantum filter ( fTSl ) forgets its 
initial condition pg' almost surely and thus admits some 
strong contraction properties as indicated by Theorem [3] 

With the truncation to n'"'" photons, convergence is proved 
only in the finite dimensional case. But feedback (|7| and 
quantum filter ( fTSl ) are still valid for = +oo. We 
conjecture that Theorems [T] and |2] remain valid in this case. 

In the experimental results reported in [10], [5], [4] the 
time-interval corresponding to a sampling step is around 
100/is. Thus it is possible to implement, on a digital com- 
puter and in real-time, the Lyapunov feedback-law (|7]) where 
p is given by the quantum filter ( fTsl l. 



Detection (-1=g, +1=e). Total number of wrong clicks (red): 10 
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Fig. 3. A single closed-loop quantum trajectory with a false detection 
probability of 1/10. 



a: P(n<3). Black: P(n=3|. Red: P(n>3) 



a: P(n<3). Black: P(n=3|. Red: P(n>3) 





Fig. 4. Averages of 10^ closed-loop quantum trajectories similar to the 
one of figure |2] (left, rjf = 0) and |3] (right, r/f = j^). 



VII. Appendix: stability theory for stochastic 

PROCESSES 

We recall here Doob's inequality and Kushner's invariance 
theorem. For detailed discussions and proofs we refer to [8] 
(Sections 8.4 and 8.5). 

Theorem 4 (Doob's Inequality): Let be a Markov 

chain on state space S. Suppose that there is a non-negative 
function V{x) satisfying ¥.[V(Xi) \ = x) - V[x) = 
— k{x), where k{x) > on the set {s : V{x) < A} = Q\. 

Then P ( sup V{Xn) > X \ Xq ^ x] < Further- 

\oo>ri>0 / 

more, there is some random t; > 0, so that for paths never 
leaving Q\, V{Xn) v > almost surely. 
For the statement of the second Theorem, we need to use 
the language of probability measures rather than the random 
process. Therefore, we deal with the space Ai of probability 
measures on the state space S. Let jio — ip he the initial 
probability distribution (everywhere through this paper we 
have dealt with the case where fiQ is a dirac on a state pq of 
the state space of density matrices). Then, the probability 
distribution of X„, given initial distribution ip, is to be 
denoted by Unif)- Note that for to > 0, the Markov property 

implies: Hn+mi^p) = Mn(Mm(<y5)). 

Theorem 5 (Kushner's invariance Theorem): Consider 
the same assumptions as that of the Theorem |4] Let pq ^ ip 
be concentrated on a state x^ E Q\ {Q\ being defined as in 
Theorem nil, i.e. (p(a;o) — 1. Assume that < k{Xn) 
in Q\ implies that Xn {x \ k{x) = Q} r\ Q\ = K\. 
Under the conditions of Theorem |4] for trajectories never 
leaving Qx, Xn converges to K\ almost surely. Also, the 
associated conditioned probability measures pn tend to the 
largest invariant set of measures M whose support set is 
in Kx. Finally, for the trajectories never leaving Qx, Xn 
converges, in probability, to the support set of M. 
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